Dipolar condensates confined in a toroidal trap: ground state and vortices 
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We study a Bose-Einstein condensate of ^^Cr atoms confined in a toroidal trap with a variable 
strength of s-wave contact interactions. We analyze the effects of the anisotropic nature of the 
dipolar interaction by considering the magnetization axis to be perpendicular to the trap symmetry 
axis. In the absence of a central repulsive barrier, when the trap is purely harmonic, the effect 
of reducing the scattering length is a tuning of the geometry of the system: from a pancake- 
shaped condensate when it is large, to a cigar-shaped condensate for small scattering lengths. For 
a condensate in a toroidal trap, the interaction in combination with the central repulsive Gaussian 
barrier produces an azimuthal dependence of the particle density for a fixed radial distance. We 
find that along the magnetization direction the density decreases as the scattering length is reduced 
but presents two symmetric density peaks in the perpendicular axis. For even lower values of 
the scattering length we observe that the system undergoes a dipolar-induced symmetry breaking 
phenomenon. The whole density becomes concentrated in one of the peaks, resembling an origin- 
displaced cigar-shaped condensate. In this context we also analyze stationary vortex states and their 
associated velocity field, finding that this latter also shows a strong azimuthal dependence for small 
scattering lengths. The expectation value of the angular momentum along the z direction provides 
a qualitative measure of the difference between the velocity in the different density peaks. 

PACS numbers: 03.75.Lm, 03.75.Hh, 03.75.Kk 



I. INTRODUCTION 

Due to the large magnetic dipole moment that 
chromium atoms possess, a gas of ultracold ^^Cr is a very 
suitable system for studying the effects of the dipolar 
interaction in different types of condensates. Although 
chromium condensates show both dipolar as well as s- 
wave contact interactions, a possible way to enhance the 
effects of the former is to reduce the value of the scat- 
tering length, which can be achieved with Feshbach res- 
onances [l|]. The wide range of configurations attainable 
in dipolar condensates have fostered an intense theoreti- 
cal and experimental activity in the field (see Ref. for 
a recent review). There have been studies delving into 
the anisotropic character of the interactions 0, 11 , the 
sensitivity of the system to the trap aspect ratio a, Q , 
the Thomas-Fermi limit in dipolar condensates [? , the 
effects of long-range in the collective excitation spectrum 
[^-[llj , the appearance of a roton minimu m fl^ lisj , and 
the instability and collapse of the system p4i - [l6| , among 
others. 

The experimental realization of dipolar Bose-Einstein 
condensates (BECs) in harmonic traps [13, [13 opens the 
possibility of studying their properties in more complex 
confining potentials. Optical lattices have been already 
addressed in the literature (see the review [2| and refer- 
ences therein), as well as double- well potentials [l9l. [20j. 
but the present work is to our knowledge the first where 
dipolar condensates are studied in toroidal traps. Exper- 



imentally, toroidal traps can be realized by shining the 
condensate with a blue-detuned laser beam [2l| . 

Spontaneous symmetry breaking in quantum gases has 
recently attracted a lot of interest, especially in double- 
well confining potentials. For contact interaction conden- 
sates, the authors in Ref. (2^ have predicted that this 
phenomenon can be produced by an external inhomoge- 
neous modulation of the scattering length. Within the 
framework of dipolar condensates, symmetry breaking 
phenomena have also been predicted to appear [l^, . 
In Ref. [l^, the authors have observed that depending 
on the magnetization direction the system undergoes a 
symmetry breaking, while in Ref. [2^ it is the strength of 
the dipolar interaction that drives the same kind of phe- 
nomenon. In all these cases the symmetry broken config- 
uration is related to the self-trapped state of a Josephson 
junction. 

In this work we are interested in determining the 
physical consequences of the combination of long range 
anisotropic interactions with a ring-shaped condensate. 
We concentrate in the particular situation in which the 
magnetization direction is perpendicular to the trap sym- 
metry axis. Due to the anisotropic character of the dipo- 
lar interaction, the particle density does not conserve 
the azimuthal symmetry of the confining potential, but 
presents two well defined peaks in the direction perpen- 
dicular to the magnetization. We find that when the 
scattering length is reduced, starting from the chromium 
natural one, the height of these peaks increases. More- 
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over, for a small enough value one of them disappears, 
giving rise to a symmetry breaking phenomenon. This 
behavior could be experimentally confirmed since either 
trapping atoms in toroidal traps [23| . controlling the s- 
wave scattering length @ and the condensation of dipolar 
gases [13, Ell ^'I'S nowadays successfully achieved. 

Another important issue that has always attracted 
much interest in the field is the superfluid character of 
BECs. The presence of quantized vortices is a clear sig- 
nature of superfluidity in quantum systems [23 | and this 
subject has received intensive research in the last decades 
(see [2^ for a recent review) . It has been recently shown 
that purely s-wave condensates in toroidal traps are also 
capable of sustaining (metastable) vortex states, which 
produce observable persistent flows [2l[. This fact has 
attracted renewed interest in the physics arising in multi- 
ply connected geometries. In particular, in the lastyears, 
the formati on l26ll . stability jlSl, dynamics 11^ H^, and 
dissipation [30| of vortices in toroidal traps have been 
investigated. 

The particular characteristics of dipolar interactions 
in BECs have raised numerous questions regarding the 
physics of vortices in such systems. Strong theoretical 
effort has been invested into the characterization of vor- 
tices in dipolar condensates [3ll - [35j |. In this work we 
study stationary vortex states whose vorticity is along 
the symmetry axis of the toroidal trap. Since the dipolar 
interaction introduces a density angular dependence, the 
velocity field of the vortices is accordingly modified. To 
gain insight into this effect, we study the velocity fields 
created by phase-imprinted vortex states as a function of 
the scattering length, i.e. of the angular anisotropy of 
the condensate density. 

This work is organized as follows. In Sec.|ll]we describe 
the theoretical framework and the system under study. 
In Sec, mil we investigate the effects on the ground state 
that appear from the combination of the dipolar interac- 
tion and the contact interaction when the trap geometry 
is changed. We work with either a harmonic or a toroidal 
trap and discuss the differences between both cases. We 
tune the scattering length to analyze the case in which 
it is small enough to be comparable to the dipole-dipole 
interaction, and show that in the toroidal trap a symme- 
try breaking phenomenon arises. In Sec. IIVI we analyze 
vortex states and their associated velocity field. Finally, 
a summary and concluding remarks are offered in Scc.lVl 



II. THEORETICAL FRAMEWORK 

We consider a Bose-Einstein condensate with N = 
5 X 10^ atoms of ^^Cr, at zero temperature, confined in 
a toroidal trap. The trapping potential Vtrap(r) is the 
sum of an axially symmetric harmonic potential and the 
potential provided by a Gaussian laser beam 



potential in cylindrical coordinates reads. 



21|. The 



laser intensity is proportional to Vq and it propagates 
along the z direction with waist wq- Thus the trapping 



l^rap(r) 



{iulr^ + uP-.z") + cxp(-2 rV wl) , (1) 



where = -I- m is the atomic mass, and cjj^ 
and W2 are the radial and axial angular trap frequen- 
cies, respectively. These harmonic trap frequencies are 
fixed at cjj^ = 8.4 x 27r s~^ and = 92.5 x 27r s~^, 
giving an aspect ratio of A = ujz/uji_ = 11 and thus de- 
scribing a pancake-shaped configuration, being z the trap 
symmetry axis. The laser parameters we have used are 
Vb = Ibhuj^ and wo = lO/im. 

The magnetic dipole moment of ^^Cr is /i = 6/is, 
where hb is the Bohr magneton. We consider that 
chromium atoms interact not only via s-wave contact in- 
teractions but also via the dipole-dipole interaction 



Wdip(r - r') 



jUoM^ 1 ^ 3 cos^ 



47r 



r — r' 
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where /xq is the vacuum permeability, r — r' is the dis- 
tance between the dipoles, and 9 is the angle between 
the vector r — r' and the dipole axis. In this work we 
choose the magnetization axis to be the y-axis, which is 
perpendicular to the symmetry axis of the trap . Since 
the dipoles are situated head to tail along y direction, the 
dipole-dipole interaction is attractive along this direction 
and repulsive in the perpendicular ones. 

For weakly interacting dipolar BECs the Gross- 
Pitaevskii (GP) equation 



2m 



trap 



-g|V(r)|2 + l-dip(r) 



/iV(r) 



(3) 



has proven to provide a good description of stationary 
states. In this equation iP{t) is the condensate order 
parameter normalized to the total number of particles 
and ji is the chemical potential. The contact interac- 
tion potential is characterized by the coupling constant 
g = Anh^a/m, where a is the s-wave scattering length. 
The mean- field dipolar interaction Vdip(r) is given by 



Fdip(r) = / dr'i>dip(r-r')|V(r')| 



(4) 



which can be more easily computed by means of 
Fourier transform techniques (see |35| and references 
therein). Inserting expression ^ into Eq. ([3]) an integro- 
differential GP equation is obtained, which we solve using 
the imaginary-time method. 

The energy density functional has the standard form 
but with the additional term E'dip: 



kin 



2m 



trap 



Edip — 

14rap|V'Prfr- 



Kiipl^I'dr. 



(5) 



The last term E^ip represents the energy due to the dipo- 
lar interaction. 
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III. GROUND STATE 

In this section we analyze the ground state of the dipo- 
lar system, ipo{r), which is found by minimizing the total 
energy ([5]), that is, by solving the GP equation ([3]). We 
separately address the cases of harmonic and harmonic 
plus Gaussian potential, to better analyze the modifica- 
tions the latter introduces. 



A. Harmonic trap 

First, we analyze the effect of the dipolar interaction 
in a BEG confined in a harmonic trap, i.e. without any 
Gaussian barrier (Vq = in Eq.JT])), and with the dipoles 
aligned along the y-axis. Since in this direction the dipo- 
lar potential is attractive, it is energetically favorable for 
the system to locate more atoms along the y-axis rather 
than in the x one. Therefore, the net size of the con- 
densate along the magnetization direction increases with 
respect to the purely s-wave case (for a more detailed 
discussion see Ref. [Sa] and references therein), whereas 
along X and z directions it decreases. When the dipolar 
effects arc magnified, that is the scattering length is re- 
duced, the anisotropic nature of the interaction becomes 
even clearer. 

From Fig.[T]it can be seen that with a scattering length 
a — 100 as the condensate is still pancake shaped, with 
the a; and y sizes almost equal and the z one much 
smaller. This shows that for large scattering lengths 
the dipolar effects are almost negligible, since the dipole- 
dipole interaction introduces only a slight perturbation 
of the density and therefore the sizes in x and y direc- 
tions are similar. When the scattering length is reduced 
the dipolar effects become enhanced and the magneti- 
zation direction is privileged, yielding therefore a larger 
size than in the perpendicular direction. For a further 
reduction of the scattering length down to a = 12 as the 
condensate resembles a cigar with a similar size in the 
X and z directions but a much larger one along y. For 
a < 12aB the system becomes unstable. 

In order to quantitatively analyze the deformation of 
the condensate in the pancake trap, we compute the root- 
mean-square radius in each direction, namely Ri where 
i = x,y^ z. Thus for a ~ 100 as we find Rx = 10.41 /xm. 



fl=100£I„ 



Ry = 10.75 /im and Rz = 1.26 ^m and thus R^ < Rx — 
Ry while for a ~ 12 as we obtain R^ = 1.36 /im, Ry = 
6.13 iim and Rz = 0.87 //m which verifies Rx — Rz < Ry 
Therefore, in a BEG confined in a pancake trap potential 
with the dipoles aligned perpendicularly to the trap axis, 
a change in the scattering length may be translated to 
a significant change in the condensate geometry. In this 
particular case, just by reducing the scattering length 
from 100 as to 12 as, it is possible to tune the geometry 
from a nearly pancake-shaped to a nearly cigar-shaped 
condensate. 

To understand this behavior, we display in Fig. [2] the 
dipolar potential given by Eq. @ for a — 20 ub- One 
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FIG. 1: (Color online) Ground state densities along a::-axis 
(solid line), j/-axis (dashed line) and z-axis (dotted line) for a 
harmonically trapped condensate for three different scattering 
lengths: a = 100, 20, 12 as (from left to right). 
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FIG. 2: Dipolar potential, Eq. Q, as a function of x and y at 
z = for a harmonically trapped condensate with a = 20 a_B . 



can see that there exists a negative minimum at the cen- 
ter of the trap and that the potential goes to zero more 
rapidly in the x direction than in the y direction. More- 
over there are two maxima around |a:| = 16.5 ^m with 
Vmax = lfi(^±- Thus the energy cost to accommodate 
atoms along the y axis is lower, which produces the con- 
densate elongation in this direction. 



B. Toroidal trap 

Now we turn to the study of dipolar condensates in 
toroidal traps. The Gaussian potential in Eq. ([T]) in- 
troduces a repulsive barrier along the z axis. For large 
scattering lengths this barrier is not enough to burn a 
hole at the center of the condensate. However, when 
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FIG. 3: (Color online) Density profiles of the ground state 
configuration along 2;-axis (solid line), y-axis (dashed line) 
and 2-axis (dotted line), for a — 100,20, 14 as (from left to 
right). 
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the scattering length is reduced the chemical potential 
of the condensate becomes smaller and eventually equals 
Vq, under which condition a hole appears at the center. 
By changing Vg one can choose the scattering length for 
the onset of toroidal geometry. In what follows we will 
work with Vq = 15 huj±, for which a hole appears slightly 
above a = 30 (see Table H]). We have checked that 
no qualitatively different behavior appears when other 
intensities of the Gaussian beam are used. 



TABLE I: Values of the root-mean-square values and chemical 
potential for different scattering lengths for Vo — 15hLjj±. 
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30 


9.65 
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20 


9.37 
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14 
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11.749 


12 (SB) 


10.24 


6.78 


0.89 


8.693 



When the scattering length is large, the main effect 
of the Gaussian barrier is to decrease the central density 
and expand the cloud, as shown in the left panel of Fig. [3] 
for a ~ 100 as. In this ease, the radii fulfill — Ry, 
with Rx slightly smaller than Ry (see Table H]). This 
behavior is similar to the previous harmonically trapped 
system because the dipole-dipole interaction introduces 
only a slight perturbation of the density. 

When the scattering length is small, however, the sit- 
uation is completely different from the purely harmonic 
trap. In this latter, the mean-field dipolar interaction was 
attractive along the y direction and repulsive along the 
X direction, resulting in a broadening of the condensate 
in the y axis, namely Rx < Ry In contrast, in a toroidal 
trap and in particular below a = 20 as, the rms-radii ful- 
fill Rx > Ry (see Table U) , which certainly presents the 
opposite behavior to the purely harmonic trap system. 

To understand the difference between the two geome- 
tries we have drawn in Fig. 2] the mean-field dipolar po- 



FIG. 4: Dipolar potential, Eq. Q, as a function of x and y 
at the z = plane for a toroidal condensate with a = 20 as. 



tential for a condensate with a — 20 as (compare with 
its corresponding harmonically trapped system in Fig.[2|). 
We can see that it presents two deep minima on the x 
axis and two saddle-points on the y axis. The main con- 
sequence of this energy landscape is that the system will 
present attractive regions in the x axis and repulsive re- 
gions in the y axis, in contrast to the harmonic trap con- 
figuration. This situation can be easily understood with 
the help of Fig. [5j In this figure we have schematically 
delimited two types of regions A and B with dashed lines. 
Within region A the dipoles are mainly displayed side-by- 
side, giving a net repulsive interaction. In region B the 
atoms are mainly located head-to-tail, which gives an ef- 
fective attractive interaction. A straightforward effect of 
the location of regions A and B is that it is energetically 
favorable for the system to accommodate more dipoles in 
regions B. Therefore, a dipolar condensate confined in a 
toroidal trap shows large density maxima in the perpen- 
dicular direction to the magnetization axis, and small 
maxima (or saddle-points) in the magnetization direc- 
tion. This is clearly seen in Fig. |31 where in the central 
and right panels the density profiles along x, y and z are 
shown for condensates with a = 20 as and a = 14 as, 
respectively. Note that in these cases the density along 
the z axis (dotted line) is zero, since the Gaussian poten- 
tial creates a hole along this axis. With these profiles in 
mind it is important to point out that the differences in 
the rms-radii (Table |T]) are mainly provided by the dif- 
ferences in the height of the density peaks and not by a 
major change in the size of the condensate itself, which 
remains around the same value (~ 20 /xm) in both direc- 
tions. 

To quantitatively analyze the azimuthal density depen- 
dence, we show in Fig. [H] the evolution of the maximum 
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FIG. 5: (Color online) Diagram qualitatively showing the ef- 
fective dipolar interaction in the torus. 
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FIG. 6: (Color online) Density values at the critical points 
{xm,0,0) and (0, ym,0), which correspond to a maximum 
(solid line) and a saddle point (dashed line), respectively. In 
the inset the quotient of these heights is drawn. 



density values along the x-axis, po{xm), and along the 
y-axis, poiVm), as a function of the scattering length. A 
strong deviation of both quantities occurs when the scat- 
tering length is reduced below the value in which the 
condensate geometry becomes multiply connected. The 
vertical line in Fig. [6] marks the scattering length below 
which a hole in the condensate density appears. 

For an even smaller value of the scattering length, be- 
low a ~ 13 as, a symmetry breaking (SB) occurs. The 
atoms are accumulated only around one of the peaks. 
We have checked this result by starting the minimization 
procedure with different initial wave functions, including 
a symmetric wave function and a randomly generated 
one. In all cases we have obtained the same final asym- 
metric state. In Fig. [7] we compare the equidensity lines 
in the z = plane of two different configurations with 
a — 20 as and a ~ 12 gb- For a = 20 qb the system 
still presents reflection symmetry but for a = 12 ub the 
symmetry is broken. It is interesting to observe that not 



FIG. 7: (Color online) Density contours for the ground state 
configuration in the (a;, y) plane for a condensate with a = 
20 as (dotted curves) and a = 12 os (solid curves). 



only do the atoms accumulate at one side of the con- 
densate but also their density distribution appears more 
resistive to curve. This effect is a consequence of the fact 
that the dipolar interaction forces the particles to locate 
themselves hcad-to-tail. 

The symmetry breaking phenomenon presented here 
is a result of the combination of the toroidal geometry 
of the external potential and the anisotropic character 
of the mean-field dipolar interaction, which yields a self- 
induced energy barrier in the azimuthal direction (see 
Fig. 2]). When the scattering length is small, the attrac- 
tive part of the dipolar interaction becomes large enough 
to be energetically favorable for the system to distribute 
all atoms in only one of the minima of the potential, 
which yields a symmetry broken configuration. With the 
present choice of parameters, a further reduction of the 
scattering length leads to the system collapse. 

Recently, SB phenomena in dipolar condensates con- 
fined in double- well traps have been reported in Refs. [l^ 
and [1^ . In Ref. [l^ , the rnagnetization direction is used 
to induce SB, while in Ref. [20] it is the number of dipoles 
in the double well that drives it. The mechanism trig- 
gering the effect is the presence of a repulsive barrier in 
systems that arc dominated by attractive interactions. 
The main difference between the present configuration 
and the ones exposed in Refs. [l^ US] is that we deal 
with a self- induced SB, since it is the dipolar interaction 
itself and its anisotropic character what brings about the 
effect. In contrast, in the double-well case, SB should be 
observed in purely s-wave condensates with attractive in- 
teractions, provided the number of atoms is smaller than 
the critical one. 
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IV. VORTEX STATES 

Toroidal traps have been experimentally shown to be 
capable of sustaining persistent flows [2lj . This fact is re- 
lated to the existence of metastable vortices [l^l • In this 
section we study and characterize the anisotropic velocity 
flow that a phase-imprinted metastable vortex introduces 
in the toroidal condensate in presence of dipolar interac- 
tions. 



A. Numerical generation of vorticity and vortex 
density profiles 

To generate vortex states along the z axis we have used 
an ansatz already described in Ref. [s^, namely we first 
obtain the ground state order parameter ■0o from Eq. ([3]) 
and then imprint a velocity field in the following form: 
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FIG. 8: (Color online) Density profiles for the ground state 
(solid line) and two-vortex state (dashed line) configurations 
along X-axis, y-axis, and z-axis (left, middle, and central pan- 
els, respectively) for a condensate with a — 100 as. 



ip{r) = ■0o(r) 



X + iy 



(6) 



with n the quantum number related to the velocity cir- 
culation. We use this ansatz as the initial wave function 
of the imaginary-time evolution of the GP equation ([3]). 
After convergence, we obtain the n-vortex state of the 
system. Note that with our choice of magnetization axis 
(y) the axial symmetry is removed and thus the angular 
momentum along the z axis is no longer = nhN, since 
the angular momentum operator does not commute with 
the Hamiltonian anymore. 

A signal that the system can sustain metastable vor- 
tices [23| is the presence of a valley in the ground state 
density landscape, which may produce a local vortex 
energy minimum. In this case, after the minimization 
process, one can obtain a vortex state captured in the 
toroidal trap due to a vortex energy barrier produced 
by the surrounding density maxima [27[. In our system, 
where the axial symmetry is not achieved, the energy bar- 
rier [231 is related to the height of the saddle points of 
the density lansdcape which are located along the j/-axis. 

In Fig. [8] we show the ground and two- vortex state den- 
sities for a dipolar condensate with a = 100 a^. Since 
the ground state is a simply connected condensate, the 
density does not exhibit a hole, hence we find that mul- 
tiply quantized vortices are less favored than multiple 
singly-quantized vortices. Although the initial imprint- 
ing method produces a doubly-quantized vortex along 
the z-axis at a; = and y = 0, during the minimiza- 
tion process the vortex splits into two singly-quantized 
vortices. Since the density maximum (saddle-point) is 
smaller along the y direction the final state converges to 
two vortices placed at the y-axis (see central panel of 
Fig. Ill) . This stationary configuration is possible because 
the velocity field due to the other vortex cancels with 
the contribution to the vortex velocity provided by the 
density inhomogeneity [STj . 
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FIG. 9: (Color online) Density profiles for the ground state 
(solid line) and a single vortex state (dashed line) configura- 
tions along a;-axis and y-axis (left and right panels, respec- 
tively) for a condensate with a = 20 as. 



We have found that the initial doubly-quantized vor- 
tex always splits in two singly-quantized vortices, even in 
the case of multiply connected geometry, i.e. with zero 
density at the center of the trap. However, for scattering 
lengths between a = 25 as and a = SOob the density at 
the center of the trap is very small, since the vortices are 
very close together. Below a = 25 qb the system is not 
capable of sustaining two vortices, since the height at the 
saddle-points is not large enough. 

To characterize the velocity field in the case of large 
asymmetric density configurations in the toroidal trap, 
we have studied singly quantized vortices. The velocity 
field will be analyzed in the next paragraph, but here 
we comment on the density distribution of such a con- 
figuration. In Fig. [5] we show the density profiles for the 
ground state and a single-vortex state for a condensate 
with a scattering length a = 20 as. Note that along the 
z-axis both densities vanish. In this case, in contrast to 
the case of a large scattering length, the presence of the 
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vortex docs not produce a large change in the density 
distribution compared to the vortex-free configuration. 



B. Vortex velocity field 

The initially imprinted velocity field, Eq. ([5]), of a sin- 
gle vortex corresponds to a homogeneous medium or an 
axially symmetric condensate. Our system is not axially 
symmetric and thus the modulus of velocity field, for a 
given radius, varies around the torus. This is simply un- 
derstood by reminding that in stationary conditions the 
current intensity is constant along the torus. 
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FIG. 10: (Color online) Velocity field components Vy (solid 
line) and Vx (dashed line) at the density maxima Xm and 
respectively, as a function of the scattering length. In the 
inset, the ratio between the two is shown. 



where dS is the surface element. The velocity field v(r) 
is related to the phase S{r) of the wave function through 



v(r) = — VS'(r) 

m 



(8) 



The integral in Eq. ([7]) is calculated in a torus section. It 
is easy to see from this equation that for the angles where 
the density is lower (regions A) the velocity is larger, 
while in the angles where the density shows a maximum 
(regions B) the velocity is minimum. 

In Fig.[TU]we display the x and y components of the ve- 
locity field at the positions j/^ and Xm of the density max- 
ima, respectively Vx{ym) and Vy{xm) (see also Fig. [TT|). 
For a = 100 as, since the dipolar effects are small the 
densities along each axis are similar and, therefore, both 
components of the velocity are equal. When we reduce 
the value of the scattering length, dipolar effects become 
sizeable and the density difference increases (see Fig. [6]) , 
hence decreasing the difference in both components of 
the velocity. This can be used to generate a desired lo- 
cal velocity field intensity by only tuning the scattering 
length using a Feshbach resonance. 

Although the angular momentum is not conserved, 
its expectation value (L^) can provide a qualitative es- 
timation of the difference between the velocity field in 
the different regions of the torus. As expected, for large 
scattering lengths the angular momentum per particle 
is almost one, see Fig. [T^ As a decreases the angular 
momentum decreases, being the variation much stronger 
when the condensate exhibits a hole, below a — 30 a^. 
The reduction in the angular momentum is an evidence of 
the presence of self-induced energy barriers, which dimin- 
ish the net particle flow. Since in our case the azimuthal 
energy barrier is larger when the scattering length is re- 
duced, a measure of the angular momentum in the z di- 
rection could provide information of the strength of the 
dipolar interaction with respect to the contact interac- 
tion. 
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FIG. 11: (Color online) Left panel: schematic diagram repre- 
senting the components of the current from which the veloci- 
ties depicted in Fig. [10] are taken. Right panel: velocity field 
on the plane z = Q for a = 20 as . 



V. SUMMARY AND CONCLUDING REMARKS 

In this paper we have studied the effects induced by the 
anisotropic nature of the dipolar interaction in a Bose- 
Einstein condensate of chromium atoms, especially in the 
case of toroidal trapping. Wc have considered the mag- 
netization axis to be perpendicular to the harmonic trap 
symmetry axis and we have analyzed the system for dif- 
ferent scattering lengths. 

We have studied the system confined in a harmonic 
potential. In this case it is possible to tune its geometry 
just by varying the scattering length: for a large scatter- 
ing length, the system shows mainly an oblate (pancake- 
shaped) geometry, while for small scattering lengths it 
can reach the configuration of a prolate (cigar-shaped) 
condensate. 

The anisotropic character of the dipolar interaction can 
be further enhanced if a toroidal trap is used. The den- 
sity then becomes inhomogeneously distributed along the 
torus, showing large and small density peaks as the az- 
imuthal angle is varied for a fixed radius. The larger 
density peaks arc observed in the x-axis, where the mi- 
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FIG. 12: Angular momentum as a lunction ol the scattering 
length. 

croscopic dipole-dipole interaction is repulsive but the 
net mean-field dipolar interaction is mainly attractive. 

The difference between these density peaks in the x 
and y axes increases as the scattering length is reduced. 
For a low enough value of the scattering length, which in 
our configuration lies below a ~ 13 a^, a dipolar-induced 
symmetry breaking phenomenon occurs: the atoms are 
concentrated along only one of the major density peaks, 
reminding thus of a cigar-shaped condensate parallel to 
the magnetization direction. 

Since toroidal traps can sustain persistent flows we 
have numerically obtained mctastable vortex states using 
a minimization procedure with an initial phase- imprinted 
wave function. For the geometry of the systems studied 



in the present work, the velocity fields show a strong 
azimuthal dependence. Following a given constant ra- 
dius along the torus, the absolute value of the velocity 
is smaller (larger) where the density shows a maximum 
(minimum), in the intersection with the x axis {y axis). 
We have also computed the expectation value of the an- 
gular momentum along the trap symmetry axis. As we 
have chosen the dipolar interaction direction perpendic- 
ular to this axis, the z component of the angular momen- 
tum is not quantized. Its reduction with respect to the 
non-dipolar system is an evidence of the presence of self- 
induced energy barriers, which diminish the net particle 
flow. Since the current intensity is conserved along the 
torus, the net consequence of this reduction is a larger dif- 
ference in magnitude of the velocity at the density critical 
points. 

In conclusion we have found that the combination 
of dipolar interactions and toroidal trapping potentials 
leads to sizeable changes in the density and velocity field 
distributions that could be in principle experimentally 
observed. 
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